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ABSTRACT 

Observations of recurrent explosive events (EEs) with time scale of 3-5 minutes are reported. These EEs 
have been observed with the Interface Region Imaging Spectrograph (IRIS) and have a spatial dimension of 
~ 1.5" along the slit. The spectral line profiles of C II 1335/1336 A and Si IV 1394/1403 A become highly 
broadened both in red as well as blue wings. Several absorption lines on top of the broadened profiles were 
identified. In addition, emission lines corresponding to neutral lines such as Cl I 1351.66 A, C I 1354.29 A, 
and C I 1355.84 A were identified. The C I 1354.29 A, and C I 1355.84 A lines were found only during the 
EEs whereas Cl I 1351.66 A broadens during the EEs. The estimated lower limit on electron number density 
obtained using the line ratios of Si IV and O IV is about 10 13 5 cm -3 , suggesting that the observed events are 
most likely occurring at heights corresponding to lower chromosphere. To the best of our knowledge, for the 
first time we have detected short-period variability (30 s and 60-90 s) within the EE bursts. Observations of 
photospheric magnetic field underneath EEs indicate that negative polarity field emerges in the neighbourhood 
of oppositely directed positive fields which undergo repetitive reconnection (magnetic flux cancellation) events. 
The dynamic changes observed in AIA 1700 A, 1600 A, C II 1330 A and Si IV 1400 A intensity images 
corresponded very well with the emergence and cancellation of photospheric magnetic field (negative polarity) 
on the time scale of 3-5 min. The observations reported here suggests that these EEs are formed due to 
magnetic reconnection and are occurring in the lower chromosphere. 

Keywords: Sun: atmosphere — Sun: transition region — Sun: chromosphere — Sun: UV radiation — line: 
profiles — Magnetic fields 


1. INTRODUCTION 

The solar atmosphere is highly dynamic changing on time- 
scales of minutes to hours. Among many, explosive events 
(EEs) are one of the prominent phenomena observed in the 
solar transition region. They were discovered by Brueck- 
ner & Bartoe (1983) using the observations recorded by the 
High-Resolution Telescope and Spectrograph (HRTS) on¬ 
board Black Brant sounding rockets. EEs are characterized 
with broad line profiles with high velocity components (~ 110 
km s -1 ), which form around 10 5 K. They have a spatial scale 
of ~ 1600 km (2 ") and life time of ~ 60 s (Dere et al. 1989). 

Using the observations recorded from the Solar Ultraviolet 
Measurements of Emitted Radiation (SUMER; Wilhelm et al. 
1995) spectrograph onboard SOHO, Innes et al. (1997b) re¬ 
ported observations of EEs in chromosphere, which revealed 
the presence of bi-directional plasma jets as predicted by the¬ 
oretical models of magnetic reconnection. It was also found 
that these events often occur in bursts lasting up to 30 min, 
whereas individual events may have typical lifetimes of about 
1-6 min (Innes et al. 1997a; Chae et al. 1998). EEs may occur 
at the same location with period around 3-5 min and may be 
triggered by waves found in the solar atmosphere (Ning et al. 
2004; Doyle et al. 2006). It has also been reported that EEs 
are preferentially located in the regions with weak and mixed 
magnetic polarity (Chae et al. 1998) and are associated with 
the canceling magnetic flux (see e.g., Muglach 2008; Huang 
etal. 2014). 

Recently, Peter et al. (2014) studied similar events using 
the observations recorded by the Interface Region Imaging 
Spectrograph (IRIS; De Pontieu et al. 2014) and found exam¬ 
ple of a bi-directional jet from the emission profiles of Si IV 
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1394/1403 A doublet. They also found absorption lines from 
cooler ions superimposed on these lines suggesting that these 
hot events are occurring in the cooler atmosphere of the Sun. 
Schmit et al. (2014) also found absorption features from a 
multitude of cool atomic and molecular lines while studying 
the broadened profiles of Si IV transition region lines dur¬ 
ing the brightening events. Huang et al. (2014) studied single 
event of EE along with underlying magnetic field evolution. 
They found the evidence of magnetic flux cancellation and 
suggested that magnetic reconnection must have taken place 
during the EE. 

In this paper, we present observations of recurrent EEs us¬ 
ing the high-resolution spectroscopic and imaging observa¬ 
tions from IRIS. We also study the evolution of underlying 
magnetic field and explore the relationship between the EEs 
and the presence of waves in the atmosphere. The paper is or¬ 
ganized as follows. We describe the observations in section 2 
and discuss the results in section 3. We have summarized the 
results and concluded in section 4. 

2. OBSERVATIONS 

Data analyzed in this work were obtained by IRIS on 2014 
March 4 between 12:39 UT to 14:37 UT using sit-and-stare 
mode. The top left and middle panels of Figure 1 shows the 
location of the IRIS slit above the slit-jaw images (SJI) ob¬ 
tained in C II 1330 A and Si IV 1400 A. The top right, bot¬ 
tom left and middle panels show corresponding images taken 
in 1600 A, 171 A and 193 A passbands of the Atmospheric 
Imaging Assembly (AIA; Lemen et al. 2012) on board the 
Solar Dynamics Observatory (SDO). The bottom right panel 
displays line-of-sight (LOS) magnetic field measurement ob¬ 
tained from the Helioseismic and Magnetic Imager (HMI; 
Schou et al. 2012), also on-board SDO. The arrows in all the 
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Figure 1 . The location of EE marked with arrows detected in slit-jaw images of IRIS 1330 A, and 1400 A, and AIA 1600 A, 171 A, and 193 A images as labeled. 
The corresponding HMI magnetic field map is also shown in bottom right panel. The vertical white continuous line on the top of IRIS slit-jaw images represents 
IRIS slit position showing that the slit is passing through the EE studied here. 


panels mark the location of the recurrent EEs. 

IRIS spectra were obtained with an exposure time of 4 s 
resulting in cadence of approximate 5 s, whereas SJI were 
obtained with an exposure time of 4 s and effective cadence 
of 15 s. In this study, we have used IRIS level-2 data pro¬ 
vided by the IRIS team. IRIS slit-jaw images from different 
filters and detectors are already co-aligned for level-2 data 1 . 
We also used the fiducial mark along the slit and SJI to correct 
any offset between them. The wobble effect due to thermal 
flexing between the guide telescope and the main IRIS tele¬ 
scope is already corrected in regular IRIS operations based on 
the orbital wobble tables (De Pontieu et al. 2014, and private 
communication, Hui Tian). This correction may still leave an 
uncertainty of about 1-2 pixels, which will not be important 
for the analysis performed over several spatial pixels. 

Data obtained from AIA and HMI have also been utilized in 
this work. IRIS and AIA observations were co-aligned using 
IRIS-SJI Si IV 1400 A and AIA 1600 A images. All the HMI 
and AIA images obtained in different filters were co-aligned 
and de-rotated with respect to AIA 1600 A image obtained at 
14:09:52 UT using the standard Solar Software (SSW) rou¬ 
tines. The obtained dataset provide an unique opportunity to 
study time evolution of recurrent EEs using both imaging and 
spectroscopic observations and to study the evolution of un¬ 
derlying magnetic field. 

3. RESULTS AND DISCUSSIONS 

Transition region EEs are identified by very broad emission 
line profiles showing non-Gaussian enhancements in both the 

1 http://iris.lmsal.com/itn26/calibration.html#coalignment-between- 
channels - and- sj i- spectra 


wings (Brueckner & Bartoe 1983). We identified one such 
small-scale bright structure (with spatial width <1.5 ") in the 
slit-jaw images of C II 1330 A and Si IV 1400A at the location 
[348.07 ", -121.44 "]. The identified EE is marked with ar¬ 
rows in all the panels in Figure 1. A corresponding animation 
is provided online, see ee.mp4. An enhancement in intensity 
corresponding to the location and time of EEs identified in 
C II 1330 A and Si IV 1400 A images is observed in AIA 
1600 A image (top right panel). The enhancement can also 
be identified in images taken in AIA 171 A and AIA 193 A 
images, though not as clearly as in AIA 1600 A. The bottom 
right panel displays magnetic field map. The arrow locates 
a small-scale (« 2" x 2") negative polarity region (average 
field strength 350 G) surrounded by the positive magnetic po¬ 
larity regions, which spatially and temporally corresponds to 
the identified EEs. The vertical white line in the IRIS slit-jaw 
images locates the IRIS slit. Fortuitously, the slit was located 
right at the location where the EE occurred allowing us to 
perform detailed spectroscopic study. Below we discuss the 
spectroscopy properties of this feature. 

3.1. Evolution of Spectral Line Profiles 

We selected a small portion of IRIS slit data correspond¬ 
ing to the spatial location the EE during the time interval 
of 14:09 UT to 14:36 UT. In Figure 2, we plot time evolu¬ 
tion of the spectral line profiles of C II 1335/1336 A doublet, 
Si IV 1394 A, and 1403 A (see panels I, J and K). An anal¬ 
ysis of the time evolution of the profiles revealed that the EE 
occurred in multiple bursts with a time period of 3-5 min and 
each burst lasted approximate for 2-3 min. Time evolution of 
event indicated by the enhancement in the intensity and width 
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Figure 2. Wavelength-time plot of the observed EE in C II 1335/1336 A doublet, Si IV 1394 A, and Si IV 1403 A spectral lines (panels I, J, and K). Typical 
spectral line profiles at various locations are shown in panels L, M, and N, where continuous lines show the line profile at the peak of event (position C in panel 
J), dot-dash line indicates profile during the quiescent time (position A in panel J), and dashed line is for intermediate time (position D in panel J). Vertical 
continuous and dashed lines in panels L, M, and N indicate most prominent and less prominent absorption lines from known ions whereas dotted lines indicate 
absorption features from unknown ions. 
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Figure 3. Emission lines from neutral atoms such as Cl I 1351.66 A (panel P), C I 1354.29 A (panel Q), and 1355.84 A (panel R) observed during EEs 
corresponding to position C of panel J in Figure 2. 
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correspond to each burst (see panels I, J and K). 

In the beginning, for the first few bursts, the enhancement in 
the intensity was by a factor of 4-12 whereas for later events 
the enhancement in intensity was about a factor of 20-40 with 
respect to the pre-EE phase at the same location (see sec¬ 
tion 3.2). Since, as can be seen in the animation 4 ee.mp4’, the 
spatial extent and locations of the bursts changed with time, 
a few bursts in the beginning and at the end could only be 
observed marginally. 

Panels L, M and N, respectively, show spectral line profiles 
of C II doublet, and two Si IV lines (1394 A and 1403 A) 
at different times. The solid lines show the spectra at the 
peak time of EE (labeled as C in panel J) and dashed line 
show the spectra for the time labeled as D in panel J. For 
comparison, we have also over-plotted the spectra at a qui¬ 
escent time (labeled as A in panel J) with dashed- dotted 
line. As can be inferred from the plots, during quiescent 
phase C II 1335/1336 A doublet line show extremely weak 
self-absorption feature (almost non-existent) in the line cen¬ 
ter, which could be due to large opacity as was pointed out 
by Peter et al. (2014). At the times when EEs have started, 
profiles of C II and Si IV have strongly broadened. This is 
essentially due to the enhancement in both the wings, which 
were sometimes observed to be asymmetric. During the bursts 
phase, self-absorption features in highly broadened C II lines 
increase sharply and bifurcate the spectral line in two (see 
panel L, solid line and dashed line). The O IV 1401 A spec¬ 
tral line show only about 50% enhancement in the intensity 
during the peak phase of the EEs, unlike the spectral lines of 
C II and Si IV, where the enhancement is much higher. 

During the initial times, when EEs were caught marginally 
by IRIS slit, profiles showed mostly blue-shifted emission (or 
blue wing enhanced). However, during the later time, profiles 
became very broad with enhancements in both red and blue 
wings. We fitted the complete C II doublet with two positive 
and two negative amplitude Gaussian functions and each Si IV 
lines with a single Gaussian function. While fitting, we used 
error bars, which were provided with IRIS level-2 data on the 
measured data numbers (DN). As EEs are believed to have 
highly non-Gaussian profiles, the reduced^ 2 and the residu¬ 
als of the fit were generally higher during the peak activity of 
the EEs. We measured Doppler velocity and width (FWHM) 
of the line with respect to the corresponding average line cen¬ 
ter position obtained at the quiescent phase (at position A la¬ 
beled in panel J). Doppler width at quiescent phase is about 
35 km s -1 for both C II and 25 km s -1 for both Si IV lines. 

Both the prominent Si IV lines show almost identical 
Doppler velocity and width variation with time. Doppler 
velocity and width measured from both the Si IV lines ex¬ 
ceeds -50 km s -1 and 150 km s -1 at the peak of activity 
around 14:25 UT. However, at the same location and time, 
C II lines measured Doppler velocity and width exceeding - 
35 km s -1 and 140 km s -1 at the peak activity time, whereas 
C II self-absorption lines measured -10 km s -1 and 65 km s -1 
respectively. Unlike the two Si IV lines, the two C II emis¬ 
sion lines do not show similar Doppler velocity and width. 
C II 1335.7 A emission line consistently shows lower Doppler 
velocity and width as compared to C II 1334.5 A. The lower 
values of Doppler shift and width in C II 1335.7 A may result 
due to the presence of strong Ni II 1335.2 A absorption line. 

From the Figure 2, we identify presence of several absorp¬ 
tion features in the emission line profiles. The most promi¬ 


nent absorption feature corresponds to Ni II 1335.20 A line 
observed on the top of broadened emission line profile of 
C II 1335.71 A. Absorption features of Ni II 1393.33 A, 
and Fe II 1403.10 A were also observed very prominently 
on the top of profiles of Si IV 1393.76 A, and 1402.77 A re¬ 
spectively. These lines are marked with vertical solid lines. 
Other less prominent absorption features were also found 
at 1334.82 A, 1335.07 A, Fe II 1335.41 A, 1336.07 A, 
Fe II 1392.81 A, Fe II 1393.21 A, Fe II 1393.59 A, 
Si IV 1393.80 A, 1394.17 A, 1402.35 A, Si IV 1402.77 A. The 
identified absorption lines are marked with vertical dashed 
lines whereas unidentified lines are marked with vertical dot¬ 
ted lines. Many of these absorption features were reported 
very recently by Schmit et al. (2014); Peter et al. (2014) and 
Yan et al. (2015). Some of them have also reported the pres¬ 
ence of self-absorption lines on the top of broadened profiles 
of C II 1334.54 A and 1335.71 A, Si IV 1393.76 A and 
1402.77 A during the explosive events. 

We have also identified several emission lines coming from 
neutral atoms such as Cl I 1351.66 A, C I 1354.29 A, 
C I 1355.84 A and O I 1355.60 A. All these lines except 
the O I 1355.60 A show presence of either self-absorption 
or absorption features (see Figure 3) during the peak activ¬ 
ity time of explosive events. To the best of our knowledge, 
this is first report of appearance of self-absorption features 
in the emission lines of neutral atoms. However, we note 
that the associated error bars on the data points are relatively 
large. Therefore, a more detailed investigation is required to 
confirm that these are indeed self absorption features. The 
C I 1354.29 A, and 1355.84 A lines were not present during 
the quiescent phase and were detected only during the EE, 
whereas Cl I 1351.66 A line was present during the quies¬ 
cent phase and broadened during the EEs. The OI 1355.6A 
was present both during the quiescent phase as well as dur¬ 
ing EEs. However, surprisingly, we didn’t find any absorption 
feature and noticeable line broadening of O I 1355.6 A line 
during any time, such as Cl I 1351.66 A, C I 1354.29 A, and 
CI1355.84A during the EEs in the IRIS spectra. 

3.2. Intensity Evolution of Explosive Events 

IRIS records both spectroscopic and imaging data of the so¬ 
lar atmosphere simultaneously. While spectra recorded with 
the slit has very limited field of view, IRIS slit-jaw imager 
(SJI) can record data in the bigger field of view. As mentioned 
in the Section 3.1, IRIS slit caught few bursts marginally. 
Thus we also use SJI images obtained in 1330 A and 1400 A 
passbands of IRIS to study the full time evolution of the whole 
EEs sequence. We also looked into the different AIA pass- 
bands to find any signatures of these events in the upper layers 
of the atmosphere. 

In top panels of Figure 4 (see movie eeJris.mpA ), we show 
an area chosen to study the intensity evolution in IRIS 1330 A 
and 1400 A passbands. Middle panel of Figure 4 displays the 
change in the intensity with time as recorded by IRIS spec¬ 
trometer. The light curves clearly reveal that the first and last 
few bursts were weak. The strongest burst were seen starting 
at around 14:24 UT. In the bottom panel of Figure 4, we plot 
intensity obtained in the chosen area with time in different 
passbands of IRIS. Recently, Martinez-Sykora et al. (2015) 
found that IRIS passbands have significant contributions from 
the continuum. However, continuum effect will be significant 
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Figure 4. Top panels: Area chosen to obtain intensity evolution in IRIS 1330 A and 1400 A passbands. Middle panel: Intensity variation with time as recorded 
by IRIS spectrometer in C II 1335/1336 A doublet, and Si IV 1403 A spectral line. Bottom panel: Intensity variation with time as obtained from IRIS 1330 A 
and 1400 A images in the chosen area. 


only when signals in the C II and Si IV lines are too weak 
to explain the presence of observed features in the 1300 and 
1400 A SJIs. In the case of EEs, signals in C II and Si IV 
profiles are very strong (stronger by a factor of 4-40), which 
are seen in the SJIs, thus, the contribution from the continuum 
may be ignored. The light curve for two IRIS passbands show 
very similar evolutionary characteristics. The light curve ob¬ 
tained from IRIS passbands indicate that total intensity of EEs 
was increasing with time in the beginning and started decreas¬ 
ing only after 14:28 UT. The strongest burst in SJI images 
were seen at around 14:28 UT, which is partially being cap¬ 
tured by the IRIS slit. It should also be noted that since the 
IRIS spectra were recorded at higher cadence as compared 
with the SJI images, we are able to see even shorter bursts 
within each EE recorded by SJI. 

To find the period of intensity variability within the EE, we 
performed wavelet analysis (Torrence & Compo 1998) on the 
intensity variation with time as recorded byClI 1335/1336A 
(see Figure 5) and Si IV 1394 A and 1403 A (see Figure 6) 
spectral lines. The top panels in both the figures show the 
variation of intensity with time starting at 14:21:28 UT. The 
bottom left panels are wavelet power spectrum (color in¬ 
verted) with 99 % confidence levels and bottom right panels 
are global wavelet power spectrum (wavelet power spectrum 
summed over time) with 99% global confidence. We used 
background (trend) subtracted intensity to obtain the wavelet 
power spectrum. Background (trend) was obtained by taking 
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Figure 5. Wavelet analysis result for the intensity variation with time (start¬ 
ing from 14:21:28 UT) as recorded by IRIS C II 1335/1336 A spectral lines. 
The top panel show the intensity variation with time. The bottom-left panel 
show the color-inverted wavelet power spectrum with 99% confidence-level 
contours, while the bottom-right panel show the global wavelet power spec¬ 
trum with 99% global confidence level drawn. The periods PI and P2 at the 
locations of the first two maxima in the global wavelet spectrum are shown 
above the global wavelet spectrum. 
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Figure 6. Same as Figure 5 but for Si IV 1394 A spectral line. 

the 10-point running average of the intensity variation. The 
power spectra obtained for both the spectral lines reveal the 
presence of short-period intensity variability in the EE light 
curves with two distinct periods of around 30 s and 60-90 s in 
addition to 3-5 min variability. We also find that even with¬ 
out background (trend) subtraction, the obtained power peaks 
at the same periods in the wavelet spectra, though with lower 
confidence level. We note here that plots shown here is cre¬ 
ated for the intensities which were obtained after taking sum¬ 
mation over the spectral profile range of the lines. However, 
we have found that the variability obtained after perform¬ 
ing wavelet analysis on Gaussian amplitudes, and product of 
Gaussian width and amplitude are essentially the same. Fur¬ 
ther, we rule out the effect of wobble, if any, on the obtained 
short term variabilities. It has been found that during a course 
of half orbit (« 48 min), IRIS pointing moves by about 17 -18 
pixels (« 3", in course of full orbit, pointing may come back 
to its original position) on the Sun (De Pontieu et al. 2014). 
This means that slit location will change by about 1 pixel in 
2.8 min, whereas, in this analysis, we found intensity variabil¬ 
ities of the order of 1 min and less within the individual EEs 
extended over several spatial pixels. Therefore, instrumental 
wobble do not have affected our analysis and henceforth the 
obtained short-period intensity variabilities. 

Repetitive nature of explosive events and jets with period¬ 
icities of 3-5 min have been reported earlier (see e.g, Ning 
et al. 2004; Doyle et al. 2006) and have been attributed to 
presence of MHD waves in the atmosphere with similar peri¬ 
ods. However, to the best of our knowledge, this is the first 
time a short-period variability with periods of 30 s and 60- 
90 s has been detected in EEs. Recently, using Hi-C data Pant 
et al. (2015) found evidence of short-period (30 s and 53-73 s) 
oscillations in braided magnetic region, providing evidence in 
favor of connection between short period waves and bursts. 
However, at this point we can not conclude whether oscilla¬ 
tions are driving burst events or vice-versa. 


3.3. Electron Density in Explosive Events 
To estimate electron density during the EEs, density sensi¬ 
tive O IV 1401 A and 1399 A line pair may be used, although 
see (Dudfk et al. 2014). Both these lines are very weak in the 
IRIS spectra. By binning over a few pixels the signal to noise 
ratio (SNR) for O IV 1401 A can be improved. However, this 


does not improve the SNR for O IV 1399 A. 

Based on the suggestions provided by Peter et al. (2014), 
we estimated electron density using O IV 1401.16 A and 
Si IV 1393.76 A (see however Doyle et al. 2013), which only 
provides a lower limit of the electron density. Since our goal 
is to get an order-of-magnitude estimate for the electron den¬ 
sity, method suggested by Peter et al. (2014) will serve the 
purpose. For this calculation, we have used the photospheric 
abundances of Grevesse & Sauval (1998) and CHIANTI ion¬ 
ization equilibrium (Landi et al. 2012). 

In order to improve the signal to noise ratio, we binned 
spectra over 7 pixel along the direction of the slit. The top 
two panels in Figure 7 show two examples of O IV 1401 A 
line profiles at two different times. As can be seen from the 
line profiles there are still not enough counts to allow to fit 
a Gaussian, we summed over the profile (wavelength range 
from 1400.94 A to 1401.50 A) and subtracted the contribu¬ 
tion of continuum to estimate the O IV 1401 A contribution. 
Although the Si IV 1393.76 A line is strong enough to fit a 
Gaussian and obtain the Gaussian integrated intensity, how¬ 
ever, for consistency we have used the same method to obtain 
the intensity as for O IV 1401 A. 

The resultant intensity variation of both the lines are plotted 
in middle panel of Figure 7. Electron number densities esti¬ 
mate during the quiescent time is around ~ 10 11,7 cm -3 which 
increases to ~ 10 13,5 cm -3 during the peak of the activity (see 
bottom panel of Figure 7). We note that since the intensi¬ 
ties derived for O IV 1401 A is the upper limit, the densities 
obtained here are essentially the lower limit. This estimated 
lower limit on the electron number density together with ap¬ 
pearance of neutral lines (C 1 1354.29 A, and C I 1355.84 A) 
during the EE strongly suggests that event probably occurred 
somewhere in the lower chromosphere, as was also pointed 
out by Peter et al. (2014) 

3.4. Comparison with AIA Observations 

To study if there were any hotter counterparts of the repet¬ 
itive EEs, we examined the different AIA passbands. Fig¬ 
ure 8 displays AIA images taken just after the peak at 14:28 
UT as seen in the bottom panel of Figure 4 in all its pass- 
bands. For complete time evolution, see movie ee.aia.mpA. 
The observations recorded using 1700 A and 1600 A pass- 
band show appearance and movement of EEs in full extent - 
marked by a rectangular box - similar to that seen in IRIS- 
SJI images, see the light curves plotted in Figure 9. The light 
curves are plotted using normalized intensities, which were 
obtained as [(I(t)-min(I(t)) = max(I(t)-min(I(t))]. However, 
the hotter channels of AIA show only a small but distinguish¬ 
able brightening that is spatially correlated to that of the EEs. 
These brightening are located in the images using a square. 
The full extent of these brightening is not clear in these hot¬ 
ter channels. Also, their is a time lag of about 5 minutes in 
the appearance of these brightening in hotter channels. This 
makes us wonder if these brightening are similar to those ob¬ 
served using SJI as well as AIA 1700 A and 1600 A images. 
We would also like to point out that coronal loops that con¬ 
nect to more southerly negative flux from the sunspot is more 
likely crossing the field of view. As pointed out by Peter et al. 
(2014), EUV radiation is strongly attenuated by overlying hy¬ 
drogen so it would not be surprising to see no signature of 
the events in AIA passbands at the same time. However, as 
time progress more energy is released, resulting in ionization 
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Figure 7. Top panel: O IV 1401 A profiles obtained at two different instants of time after binning over 7 pixels along the slit at the EE location. Middle panel: 
Intensity variation of O IV 1401 A and Si IV 1393.76 A lines with time. Bottom panel: Temporal evolution of estimated electron number density (lower limit) 
obtained from the ratios of Si IV to O IV lines. 


of hydrogen with time, we see response in AIA passbands for 
later events. 

In order to understand the nature of the small brightening 
seen in the AIA hotter channels, we performed differential 
emission measure (DEM) analysis to obtain the differential 
emission measure distribution using the regularized inversion 
method of Hannah & Kontar (2012). The DEM was obtained 
using AIA 131 A, 171 A, 193 A, 211 A, 335 A, and 94 A in¬ 
tensities recorded at nearby time of 14:29:54 UT when emis¬ 
sion in AIA 193 A channel peaks. The obtained DEM curve 
is plotted in Figure 10. The red curve is DEM before back¬ 
ground subtraction and the blue curve is after the background 
subtraction. Background intensities are intensities obtained 
at the EE location just before the brightening start appear¬ 
ing in AIA images. As expected, the overall DEM values 
decreases after the background subtraction. The DEM curve 
show rather a strong dip at ~ 1 MK after the background sub¬ 
traction. However, the peak emission is still coming from a 
temperature of log T = 6.3. This suggests that the brighten¬ 
ing seen in the AIA channels is at coronal temperature. There 
is also another peak at log T = 5.8. However, the error bars 
are much bigger at this location in the plot. 

If the brightening seen in hotter AIA channels are related 
to those seen in IRIS-SJI as well as AIA 1600 A and 1700 A 
images, then the plasma must be heated to about 2 MK. As 
mentioned earlier, the time analysis shows that there is a lag 
about 5 minutes in the appearance of these brightening in hot¬ 


ter channels. Also we note that by the time these brightening 
appear in hotter channels, they disappear from the SJI images. 
However, at this point, we are not able to conclude if these 
brightening are exactly related to each other. 

3.5. Magnetic Field Evolution of Explosive Event 

We analyzed HMI/SDO magnetogram data to study evolu¬ 
tion of photospheric magnetic field during the EE. Figure 11 
(see movie eeJimi.mpA) shows presence of small-scale neg¬ 
ative field (enclosed with a rectangular box) surrounded by 
the positive field at the location of EE. We have studied the 
temporal evolution of positive and negative magnetic field as 
well as flux in the the boxed region and plotted in Figure 11. 
As can be seen from the plots, negative polarity field and flux 
(solid lines in bottom panels) increases and decreases during 
the EEs on the time scale of 3-5 min. Whereas positive polar¬ 
ity field and flux (dashed lines in bottom panels) show some¬ 
what decay during that time. 

To find any correlation between fluctuations in magnetic 
field measured from HMI and in intensity during EE mea¬ 
sured from AIA and IRIS, we perform correlation analysis at 
different time delays. HMI magnetogram is obtained at photo- 
spheric height whereas AIA 1700 A is formed at photosphere 
and temperature minimum height, thus, are more suitable for 
direct correlation study. Correlation of intensity fluctuations 
obtained from AIA 1600 A, IRIS 1330 A, and 1400 A are per¬ 
formed with respect to that obtained from AIA 1700 A. As ca- 
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Figure 8. Images obtained from AIA 1700 A, 1600 A, 304 A, 131 A, 171 A, 193 


i - 1 - 1 -■- 1 -■- 1 - 1 -r 



14:20 14:24 14:28 14:32 

Start Time (04-Mar-14 14:16:16) 


Figure 9. Intensity light curves obtained from AIA 1700 A, and 1600 A 
passbands as labeled at the location of EEs. 

dence of HMI is 45 s o whereas that of AIA 1600 A, 1700 A and 
IRIS 1330 A, 1400 A are 24 s and 15 s respectively. We inter¬ 
polated AIA 1700 A with respect to HMI time sequence using 
spline routines. We also interpolated AIA 1600 A, and IRIS 
1330 A, 1400 A light curves with respect to time sequence 
of original AIA 1700 A. Thus, we obtained cross-correlation 


A, 211 A, 335 A, and 94 A passbands as labeled during the EE. 



io 17 P , , ,_,_i_,_,_,_,_i_,_, ,_,_i_,_ 

5.0 5.5 6.0 6.5 

log 10 T 

Figure 10. DEM curve obtained using coronal emission data from AIA 
131 A, 171 A, 193 A, 211 A, 335 A, and 94 A passbands at the spatial lo¬ 
cation of EE and at nearby time of 14:29:54 UT. Continuous line red curve 
is obtained using the original intensity whereas dashed line blue curve is ob¬ 
tained after subtracting background. 

coefficients at different time delays for the pairs of HMI mag¬ 
netogram (negative polarity) and AIA 1700 A light curve, 
AIA 1700 A and 1600 A, AIA 1700 A and IRIS 1330 A, and 
AIA 1700 A and IRIS 1400 A light curves. We plot respec- 




















































Observation of Recurrent Explosive Events 


9 


14:14:43 UT 14:19:58 UT 14:25:13 UT 14:32:43 UT 



0 5 10 15 200 5 10 15 200 


5 10 15 200 5 10 15 20 



14:12 14:16 14:20 14:24 14:28 14:32 14:36 14:12 14:16 14:20 14:24 14:28 14:32 14:36 

Stort Time (04-Mar-14 14:09:28) Start Time (04-Mar-14 14:09:28) 


Figure 11. Time evolution of underlying photospheric magnetic field of EE as recorded by HMI magnetogram. Top panels show the small region where EEs 
were observed. Contour levels are over-plotted at the field strength of 0 G. Over-plotted box indicate area chosen to estimate average field strength and flux at the 
EE source region. Bottom panels provide the time evolution of positive (dashed line) and negative (continuous line) field strength (left), and flux (right). 



Figure 12. Cross-correlation analysis between intensity evolution of EEs obtained from AIA 1700 A (continuous line, all panels) with HMI magnetogram 
(negative polarity, top panel, dashed line), AIA 1600 A (second panel, dashed line), IRIS 1330 A (third panel, dashed line), and IRIS 1400 A (bottom panel, 
dashed line). Cross-correlation co-efficient values obtained with original signal at different time delays are plotted in the right panels for the respective light curve 
pairs. 
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tive pairs of light curves and correlation coefficients at differ¬ 
ent time delays in Figure 12. Maximum correlation value be¬ 
tween HMI magnetogram (negative polarity) and AIA 1700 A 
light curve is about 0.77 at the time delay of 135 s, which 
corresponds to 3 HMI time frames. However, that between 
AIA 1700-1600 A, AIA 1700-IRIS 1330 A, AIA 1700- 
IRIS 1400 A pairs are about 0.84, 0.82, 0.79 respectively at 
the time delay of 0 s. Maximum cross-correlation co-efficient 
between AIA 1600 A and IRIS 1330 A and 1400 A are above 
0.90 at time delay of 0 s (not shown here). Thus, the obtained 
results suggest that intensity fluctuations as recorded from 
AIA 1600 A and 1700 A and IRIS passbands are connected to 
fluctuations in photospheric magnetic fields as recorded from 
HMI with the time delay of about 135 s. Increase in flux and 
held could be related to flux emergence, whereas decrease in 
magnetic held and flux could be associated with magnetic flux 
cancellation resulting due to magnetic reconnection. Mag¬ 
netic flux cancellation events underneath EEs had been previ¬ 
ously reported by Chae et al. (1998) and recently by Huang 
et al. (2014). However, to the best of our knowledge, this is 
first time we report such a correlated change in AIA 1700 A, 
1600 A, C II 1330 A, and Si IV 1400 A intensity with pho¬ 
tospheric magnetic held (in this case negative polarity) on the 
time scale of 3-5 min during the recurrent bursts of EEs. 

In order to get further insight of the magnetic held evolu¬ 
tion, we performed a long term study of the sunspot region i.e. 
from the time when it hrst emerged on the east limb on Feb 25, 
2014 AR 11990 till the time of the analyzed EEs using HMI 
observations. The region appeared as a simple sunspot (with 
negative polarity) within which an emergence of a positive 
held region was detected. This positive held region evolved 
with time and developed as a complete sunspot with positive 
polarity, suggesting a highly complex held evolution and for¬ 
mation of a delta-sunspot (Gupta et al. in preparation). With 
the evolution of the sunspot, various moving magnetic fea¬ 
tures (MMFs) were observed. Therefore, in the context, it 
is plausible to conclude that the active region is relatively a 
young with several small-scale MMFs around the sunspot. 
The EEs studied in current paper was related to one of the 
MMFs. 

4. SUMMARY AND CONCLUSIONS 

In this paper we studied observations of recurrent EEs 
using simultaneous spectroscopic and imaging observations 
recorded by IRIS, AIA images and HMI magnetogram. To 
the best of our knowledge, this is hrst report of such recurrent 
EEs using IRIS data. The recurring time scale of these EEs 
were about 3-5 min. During the event, line prohles of C II 
and Si IV showed enhanced broadening with Doppler velocity 
and width exceeding more than -50 km s -1 and 150 km s -1 re¬ 
spectively. In addition, we identihed several absorption lines 
on the top of broadened emission lines of C II and Si IV. 
Moreover, we also found a few neutral atom lines such as 
Cl 1 1351.66 A, Cl 1354.29 A, and C 1 1355.84 A with possi¬ 
ble self- absorption features. While C 1 1354.29 A, 1355.84 A 
lines appeared only during the EEs, the Cl 1 1351.66 A, which 
was present before the EEs, showed broadening during the 
EEs. The lower limit on electron densities obtained using 
the method proposed by Peter et al. (2014) was about 10 13,5 
cm -3 . Using the high cadence spectroscopic observations we 
also discovered short period variability (~ 30 s and 60-90 s) 
within the EE bursts. 

The analysis the LOS of photospheric magnetic held mea¬ 


sured by HMI underneath explosive events indicated emer¬ 
gence as well as cancellation of magnetic flux. The nega¬ 
tive polarity magnetic flux showed continuous increase with a 
periodic fluctuations, suggesting localized cancellations (see 
bottom panels of Figure 11). The changes in AIA 1700 A, 
1600 A, C II 1330 A, and Si IV 1400 A intensities correlated 
extremely well with the changes in the negative polarity mag¬ 
netic flux on the time scale of 3-5 min during the recurrent 
bursts of explosive events, which has not been reported ear¬ 
lier, to the best of our knowledge. 

The observations of strong broadening in C II and Si IV 
spectral lines along with self-absorption and estimated high 
electron density suggests that these are lower chromospheric 
features. The correlation with canceling magnetic features 
supports the idea of formation of these features due to mag¬ 
netic reconnection. We believe that these features can be ex¬ 
plained by resistive emergence of magnetic flux underneath 
the photosphere and expanding in the upper layers (see e.g., 
Pariat et al. 2004; Isobe et al. 2007). This concept has suc¬ 
cessfully explained the formation of Ellerman Bombs, which 
were originally detected in H a observations (Ellerman 1917). 
In this scenario, a number of Q loops rise due to well known 
Parker’s instability from underneath the photosphere and ex¬ 
pand. While emergence and expansion, they interact with 
each other and the process of magnetic reconnection occurs. 
The process of reconnection may produce the EEs and heat 
the plasma locally. This process may result in a cool mate¬ 
rial being stacked upon locally heated material, thus resulting 
in a cool absorption lines superimposed on the hot emission 
lines as proposed by Peter et al. (2014). Therefore, the ob¬ 
served EEs are most likely the upper atmospheric counter¬ 
part of Ellerman Bombs. This concept, however, needs to be 
throughly verified using observations recorded in H a line core 
and wings along with IRIS observations and forward model¬ 
ing of the spectral lines using MHD simulations of realistic 
solar atmosphere. 
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